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We propose to take advantage of the properties of long Josephson junctions to realize a frequency 
variable Fabry-Perot Filter that operates in the range 100 — 500 GHz with a bandwidth below 1 
GHz. In fact, we show that it is possible to exploit the tunability of the effective impedance of the 
Josephson component, that is controlled by a dc bias, to tune, up to 10% of the central frequency, the 
resonance of the system. An analysis of the linearized system indicates the range of operation and 
the main characteristic parameters Numerical simulations of the full nonlinear Josephson element 
confirm the behavior expected from the linear approximation. 
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I. INTRODUCTION 

Layered structures have been introduced in microwave 
and quasi-optical devices to produce filters with a narrow 
transmission bandwidth [l| and as filter in the high fre¬ 
quency domain 0. Layered structures, periodic or ape¬ 
riodic, have also been proposed for fundamental physics 
applications with sophisticated optics Q and microwave 
applicators 0. A frequency band can be selected in¬ 
troducing, in the terminal points of the periodic struc¬ 
tures, suitable defects; the structures thus constitute a 
sort of Fabry-Perot filter in the optical and quasioptical 
band that alternate periodically two different sections H 
and L of high and low distributed capacitance, respec¬ 
tively. Fabry-Perot like devices provide high precision 
measurements set-up for demanding applications [^, @ . 
The defect ordinarily consists in a linear Transmission 
Line (TL), thus a well known (and commonly employed) 
system can be modelled as a periodic array of A/4 TL 
that introduces identical delays, while the defect consists 
in still another TL of different length. The periodic struc¬ 
ture determines a band gap; if one introduces a defect, a 
narrow resonance appears in the bandwidth. A natural 
limit of this technology is that the resonant frequency is 
determined by the physical characteristics of the defects; 
therefore, to tune the central frequency one should build 
(and insert) a different defect. Therefore, to overcome 
this complication superconducting tunable filters (with 
lumped parameters) have been explored to realize elec¬ 
tronically controlled devices without mechanical tuning 
In this context, we propose to employ as defects 
Josephson Junctions Transmission Lines (JJTL) flol - [l^ 
to obtain a tunable filter in the sub-THz region. The 
motivation is twofold: Josephson Junctions (JJ) are very 
fast superconducting elements, suitable to build high fre¬ 
quency and low noise transmission lines [TsI - fTsjl , pseudo¬ 
cavities [l^, and generators flTl - f^ capable to perform 
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even near the THz region. Moreover, JJ’s are nonlinear 
element whose (effective) inductance can be tuned by an 
external dc bias [I2,[2lj. The latter property of JJ allows 
to control the effective inductance, at high frequency, 
with a change of the direct current through the super¬ 
conducting element (^ . [^ . In this work we examine 
how the main properties of such electronically controlled 
defect can be exploited to design a tuinable filter. The 
aim is to show that one can tune some properties of the 
device with a (relatively) simple change of the applied 
current. The main superconducting TL Josephson ele¬ 
ment consist of a so-called long JJ [iflllil. Design capa¬ 
bility for JJ are highly developed, as demonstrated by the 
recent upsurge of superconducting metamaterials 24, [2^ , 
especially for quiet quantum measurements [l^, l26j| . It is 
therefore conceivable that the proposed concept can be 
considered for practical purposes. The work is organized 
as follows. In Sect. m we retrieve, in the ABCD matrix 
formalism the dependence of the JJTL as a function 
of the dc bias. In Sect. imi we show in the linear approx¬ 
imation the main properties of the proposed structure: 
bandwidth and tunability. We also numerically investi¬ 
gate the full nonlinear JJTL and compare its behavior to 
the approximated linear analysis. Sect. IIVI concludes. 


II. MODEL 


The constitutive equations for the JJTL element of 
Fig. [T] read [^ : 
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where h is the Planck constant, e is the elementary 
charge, V and I are the distributed voltages and currents, 
respectively; Cj , Gj , and Ic are the capacitance, conduc¬ 
tance and maximum (or critical) current per unit length 
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of the superconductive element; Lq is the geometric in¬ 
ductance per unit len gth. Eqs. mm amount to a sine- 
Gordon equation [10, [ll|- The variable (j) is the gauge 
invariant phase difference between the superconducting 
wave functions that determines the Josephson super¬ 
current Ic sin(</)), the essence of the Josephson effect, as it 
corresponds to the tunneling of Cooper’s pairs. Finally, 
Eq.(|31) is the phase-voltage Josephson relation. From the 
standpoint of the electromagnetic propagation it corre¬ 
sponds to a nonlinear distributed inductance 



This is the key property of the Josephson effect that we 
want to exploit: the effective inductance Lj can be tuned 
by the external dc bias Ib, in analogy to fluxometers 
(SQUIDs) or other microwave components [^ . 
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FIG. 1. Schematic top view (not in scale) of the planar mi¬ 
crostrip circuit representing a chain of A/4 standard TL (solid 
lines rectangles) joined to JJTL defects (dashed lines rectan¬ 
gles). The ac component travels along the z-axis, for the space 
between the sections is only shown to distinguish the various 
elements. A dc current bias is only fed through the JJ ele¬ 
ment. In the following of the paper we take the number of 
TL sections Nl = 8. 

Linearizing the JJTL Eqs. CEEl) and using the dis¬ 
tributed inductance Q in the frequency domain we find: 

^ = -juiLoI (5) 

dz 

where V and / are the voltage and current phasors and 
we use the substitution rule d{ )/dt —>• joj. The fac¬ 
tor Ceq in Eq.® amounts to an equivalent condenser 
Ceq = Cj — l/{uj‘^Lj) J- Gjjjoj, that depends on the dc 
bias Ib through Eq. ([3]). Through the whole paper we use 
the term frequency as a synonymous of angular velocity. 
The linear analysis of the TL gives the (secondary) pa¬ 
rameters /3 and Z,,'. 
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In the ABCD matrix formalism the parameters 
P and Zc are fundamental to derive the input-output 
matrix relations. For instance, the parameters P and Zc 


of Eas.dTia can be employed to derive the input-output 
relations of a JJTL: 

(v,\_( cos m jZ, sin {Pl) \(V2\ 

\h )~\ jZp^ sin [Pl) cos {Pl) )\l 2 

In the following we shall indicate the matrix of the above 
Eq.® with b. 


TABLE 1. We have chosen the main physical parameters 
of the JJTL component close to well established technology 
[H, El] (we refer to Fig. [T| : I and d are the length of the 
JJTL element and of the high capacitance nonsuperconduc- 
tive TL, respectively (the length of the low capacitance TL is 
obtained by the A/4 constraint). Rj is the JJ resistance, Jc 
the Josephson critical current density, h and w the thickness 
and the width, respectively, while Cr is the relative dielec¬ 
tric constant of the insulating barrier. The Josephson physics 
[lol . [Till dictates the other parameters: IcRj = A (A is the 
superconducting gap) and Lq = goh/w. Ch and Cl are the 
capacitance per unit length of the H and L lines, respectively. 
We also denote with Lind the inductance per unit length of 
both H and L nonsuperconductive TL. The input voltage 
source amplitude is Vo = 50fiV, series connected with a gen¬ 
erator impedance Rq = 0.2914, while the final edge of the 
array is terminated by a real resistance Rl ~ 0.2914. These 
figures are but an example of a realistic conhguration with 
acceptable performances. With these parameters the usual 
normalized Josephson parameters read: coj ~ lOOGHz, the 
Josephson penetration depth is Xj ~ 40fj,m, the normalized 
dissipation a ~ 0.006. 
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The defect described by Eq. ® has the structure of an 
ordinary defect [l| that introduces an appropriated phase 
delay of the signal between the input and the output 
points [13 ■ This defect can be combined with standard 
input-output relations for A/4 elements [l^]. If L and 
H are two transmission lines with different capacitances 
(Cl and Ch for the low and high pieces, respectively) and 
with the same inductance Lind^ one obtains a band pass 
behavior due to the periodic structure. For the complete 
system (including the JJTL defect) of Figin the trans¬ 
mission matrix reads: 

T = LHLH...LHLH D LHLH...LHLH D Rl 
^ ^ ^ ^ ^ 

Nl times Nl times 

= b (LH)^^ b Rl (10) 

that completely describes the linear behavior (here, a ter¬ 
mination is described by the resistive matrix Rl) ■ If 
we define T as the result of the matrix multiplication in 
Eq.® 


T = 


A’ B' \ 
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the input impedance Zi reads: 


y, + B' 

U CRl + D^' 


( 12 ) 


Finally, denoting with Rq the impedance (purely resis¬ 
tive) of the feeding TL, the transmission coefficient F^ 
reads: 


Zj — Rg 

Zi -\- Rg 


(13) 


Equation (m, together with Eq.Q that is contained 
in the matrix D, states that the transmission coefficient 
exhibits a dip that depends upon the dc bias current R- 
This is the main result of this paper, illustrated in the 
next Section. 




Bias Current RIR 



FIG. 3. Effect of the bias current 7 = h/Ic (i.e., the bias 
current R normalized to the Josephson junction critical cur¬ 
rent 7c) on two main parameters of the Fabry-Perot filter: the 
central frequency of the passing window and the linewidth of 
the band (see Fig. [2]). The physical parameters are given in 
Table [T] with the same letter. 



Frequency [THzJ 

FIG. 2. (color online) (a) Absolute value of the reflection 
coefficient Fi as per Eq. m for a TL interrupted by a JJ 
defects (black solid line) and of a TL without defect (red 
dashed line). The frequency, or more precisely the angular 
velocity, is measured in Hertz. The bias current reads 7 = 
h/Ic = 0.8. (b) Absolute value of the reflection coefficient 
Fi for different values of the current 7 . The other physical 
parameters are given in Table U 


III. RESULTS 

To characterize the proposed Fabry-Perot filter, we in¬ 
clude an input impedance and a load to obtain relevant 
quantities as the transmission coefficient. For simplicity, 
we consider a simple voltage source and a resistive load 
at the beginning and the end of the block diagram de¬ 
scribed in Sect. m For the parameters of Table U the 


reflection coefficient IF^I as a function of the frequency is 
displayed in Fig. Oi at two values of the bias. We stress 
that the resonance dip in Fig. ^ can be tuned by a sim¬ 
ple change of the bias current 7 = Ib/R in Eq.d^]), as 
shown in Fig. [2 }d. It appears that the central frequency 
can be tuned, within the allowed range of bias current, of 
as much as 3.5GHz (i.e., about 2% of the resonance fre¬ 
quency, for this parameters choice). This is shown in Fig. 
|3l the current R can be tuned from 0 (as negative value 
do not change the behavior of the device) to a maximum 
current that is the critical current of the Josephson ele¬ 
ment (above such current the junction develops a voltage 
and the propagation of the signal is not described by the 
present model). The Figure displays the effect of the bias 
for the two sets of parameters of Table U to demonstrate 
the generality of the mechanism. 

A key point for the validity of the approach of Sect. [Ill 
and hence of the tunability described in Fig. [Sj is that 
the linearization of the inductance employed in Eq. ([4]) is 
valid. This has been checked for the values of Table |ll(a) 
in Fig. m The comparison shows that the linear analysis 
correctly captures the main effect, as it reproduces the 
input-output voltage ratio trends. The discrepancy is 
about 0.5% at low bias, to reach a maximum of about 
1.5% close to the critical current. This is to be expected, 
for the approximation ([4]) is based on a linearization of 
the sinusoidal function in Eq.Q, that best performs at 
low bias. It is important to underline that the comparison 
is performed without any fitting parameter. Also, the 
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nonlinear model is sensitive to the amplitude of the input 
signal. In Fig. 0] we have assumed an input signal of 
amplitude Vq = 0.05mV. 

The input voltage induces a current of the order 
Vo/{Rg + Zc) — O.OSmA; this value should be compared 
to the critical current of the junction IJ — 0.1mA. We 
thus expect that the linear regime breaks up for larger 
values, see Fig. [ 5 ] 
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FIG. 4. (Color online) Comparison of the input-output rela¬ 
tion, see the inset, of the JJ transmission line (red, dashed) as 
per Eq. m with the linear prediction of the matrix D obtained 
through the approximation Q (black, solid). The relative er¬ 
ror of the linear approximation is displayed as red solid. The 
input voltage reads Vb = O.OSmF. The other physical param¬ 
eters are given in Table IJa). 


IV. CONCLUSIONS 

We have shown that a superconducting metamaterial 
containing long Josephson junctions, inserted in a 
chain of A/4 layers as a discontinuity, amounts to a tun¬ 
able defect that moves the resonance characteristic fre¬ 


quency of a Fabry-Perot filter. The change of the fre¬ 
quency, as in series JJ discrete TL is achieved sweep¬ 
ing the dc current bias. The effect can be regarded as 
relevant for applications in superconducting electronics 
and devices, for it allows to tune a Fabry-Perot filter 
with an externally dc bias. Moreover, we have found 
that a linear analysis correctly captures the main effect 
(the nonlinear corrections are below 1.5%). Some words 
of caution are in order, however. In the first place, we 
have overlooked fabrication limits and parameters toler¬ 
ance, assuming that one can effectively produce the de¬ 
sired JJ. In the second place, we have fully examined the 
core part of the device, the nonlinear (superconducting) 
Josephson transmission line and the linear array, while a 
real device consists of other parts (e.g., control electron¬ 
ics to sweep the JJ bias current). Finally, the analysis 
has been performed with a single harmonic component 
(higher order expansions can be included as in Ref. El) 
and neglecting noise. The latter requires some attention, 
for the combination of noise and an ac term can induce 
undesired switches to finite dc voltage states [3l|. 



Bias Current //,//<; 

FIG. 5. The effect of the input voltage amplitude on the 
input-output relatious, see the inset. The other physical pa¬ 
rameters are given in Table [Ja). 
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